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Relativistic  quantum  mechanics  predicts  that  when  the  charge  of  a  super-heavy 
atomic  nucleus  surpasses  a  certain  threshold,  the  resulting  strong  Coulomb  field 
causes  an  unusual  “atomic  collapse”  state  which  exhibits  an  electron  wave  function 
component  that  falls  toward  the  nucleus  as  well  as  a  positron  component  that 
escapes  to  infinity.  In  graphene,  where  charge  carriers  behave  as  massless 
relativistic  particles,  it  has  been  predicted  that  highly  charged  impurities  should 
exhibit  resonances  corresponding  to  these  atomic  collapse  states.  We  have 
observed  the  formation  of  such  resonances  around  artificial  nuclei  (clusters  of 
charged  calcium  dimers)  fabricated  on  gated  graphene  devices  via  atomic 
manipulation  with  a  scanning  tunneling  microscope  (STM).  The  energy  and  spatial 
dependence  of  the  atomic  collapse  state  measured  using  STM  revealed  unexpected 
behavior  when  it  is  occupied  by  electrons. 


It  has  long  been  predicted  that  relativistic  quantum  effects  can  render 
atoms  unstable  when  the  charge  of  a  nucleus  (Z)  exceeds  a  certain  criti¬ 
cal  value  Zc  (7,  2).  In  this  “supercritical”  regime,  the  previously  stable 
atomic  bound  states  are  predicted  to  “dive”  into  the  positron  continuum 
and  acquire  a  finite  lifetime  (3).  Quite  unlike  the  standard  Bohr- state 
orbitals  seen  in  atomic  physics  and  semiconductor  impurity  systems, 
such  supercritical  quantum  eigenstates  represent  the  quantization  of  a 
semi-classical  trajectory  having  an  electron-like  component  spiraling 
inward  toward  the  nucleus  (hence  the  term  “atomic  collapse”)  before 
spiraling  back  out  and  coupling  to  a  positron-like  component  that  propa¬ 
gates  away  from  the  nucleus  (7,  3-5). 

Some  efforts  to  observe  atomic  collapse  states  have  been  made 
through  colliding  heavy-ions  and  searching  for  the  positron  created  when 
a  “collapsing”  electron  is  tom  from  the  vacuum,  but  these  results  remain 
ambiguous  due  to  the  enormously  large  Zc  ~  170  that  is  required  ( 6 ,  7). 
Hopes  to  realize  atomic  collapse  states  rose  when  theory  predicted  that 
highly  charged  impurities  in  graphene  could  display  the  same  atomic 
collapse  physics  as  atomic  nuclei,  but  with  holes  in  the  graphene  playing 
the  role  of  positrons  (8-10).  Since  charge  carriers  in  graphene  behave  as 
relativistic  particles  with  a  large  fine  stmcture  constant  value  (77-73), 
the  supercritical  charge  threshold  is  substantially  lower  than  for  atoms, 
with  Zc  on  the  order  of  unity  (8-10).  The  atomic  collapse  state  near  a 
Coulomb  impurity  in  graphene  should  appear  abmptly  as  the  impurity 
charge  is  increased  above  Zc  and  manifest  as  a  spatially  extended  elec¬ 
tronic  resonance  whose  energy  lies  just  below  the  Dirac  point.  Such 
resonances  correspond  to  the  electron-like  part  of  the  atomic  collapse 
wavefunction  (in  contrast  to  the  far  field  measurements  of  out-going 
positrons  performed  in  heavy-ion  collision  experiments)  (8-10).  Recent¬ 
ly  subcritical  Coulomb  behavior  (i.e.,  where  Z  <  Zc  and  atomic  collapse 


resonances  are  absent)  was  observed  for 
charged  impurities  in  graphene  (14), 
but  the  observation  of  atomic  collapse 
around  supercritical  impurities  has 
remained  elusive  due  to  the  difficulty  of 
producing  highly  charged  impurities. 

Here  we  report  the  observation  of 
supercritical  Coulomb  behavior  in 
atomically-fabricated  “artificial  nuclei” 
assembled  on  the  surface  of  a  gated 
graphene  device.  These  tunable  charge 
centers  were  synthesized  by  pushing 
together  ionized  calcium  (Ca)  impuri¬ 
ties  using  the  tip  of  an  STM,  thus  al¬ 
lowing  creation  of  supercritical 
Coulomb  potentials  from  subcritical 
charge  elements.  STM  spectroscopy 
was  used  to  observe  the  emergence  of 
atomic-collapse  electronic  states  ex¬ 
tending  further  than  10  nm  from  the 
center  of  artificial  nuclei  in  the  super¬ 
critical  regime  (Z  >  Zc).  By  tuning  the 
graphene  Fermi  level  (EF)  via  electro¬ 
static  gating  we  observed  behavior  in 
this  state  that  was  unanticipated  by 
theory  (8-10). 

We  transferred  chemical- vapor- 

deposition  (CVD)  grown  graphene  (75) 
onto  a  boron  nitride  (BN)  flake  placed 
on  a  Si02/Si  surface  (the  doped  Si  sub¬ 
strate  comprised  a  gate  electrode).  The 
BN  substrate  reduced  the  charge  inho¬ 
mogeneity  of  graphene  (16-18)  suffi¬ 
ciently  that  we  could  observe  the 
intrinsic  electronic  behavior  of  pristine 
graphene  adjacent  to  charged  impurities  and  artificial  nuclei.  We  used 
Ca  dimers  as  charge  building  blocks  for  constructing  supercritical  artifi¬ 
cial  nuclei  on  graphene  because  they  were  more  easily  manipulated  with 
an  STM  tip  than  Ca  monomers.  The  dimers  formed  spontaneously  by 
allowing  the  monomers  to  diffuse  on  the  surface  by  warming  to  16  ±  3  K 
for  1  to  2  min  [details  of  the  assembly  and  characterization  are  available 
as  supplementary  materials  (fig.  SI)  (79)]. 

The  effect  of  a  charged  Ca  dimer  on  nearby  Dirac  fermions  in  gra¬ 
phene  was  determined  by  measuring  both  the  electron-like  and  hole-like 
local  density  of  states  (LDOS)  near  a  single  dimer  through  STM  dlldV 
spectroscopy  performed  at  different  lateral  distances  from  the  center  of  a 
single  Ca  dimer  (Fig.  1A).  These  spectra  exhibited  a  ~130-meV  wide 
gap-like  feature  at  EF  caused  by  phonon-assisted  inelastic  tunneling,  and 
an  additional  minimum  around  Fs  =  +0.23V  associated  with  the  Dirac 
point  (20).  These  spectra  also  display  a  systematic  electron-hole  asym¬ 
metry  in  that  the  LDOS  intensity  of  states  above  the  Dirac  point  in¬ 
creased  as  the  tip  neared  the  dimer  while  the  LDOS  of  states  below  the 
Dirac  point  correspondingly  decreased.  This  asymmetric  behavior  is 
characteristic  of  subcritical  impurities  and  indicates  that  Ca  dimers  were 
positively  charged  on  graphene  at  this  doping  level  (14).  Such  electron- 
hole  asymmetry  was  observed  for  Ca  dimers  at  different  gate  voltages  in 
the  range  -60  V  <  Fg  <  +30  V  [see  supplementary  materials  fig.  S6  (79)] 
and  no  signature  of  charge  state  instability  was  seen,  indicating  that  Ca 
dimers  were  charge-stable  within  our  measurement  conditions. 

Artificial  nuclei  containing  up  to  five  Ca  dimers  were  assembled  by 
using  atomic  manipulation  techniques  to  push  individual  dimers  into 
small  clusters  (see  insets  to  Fig.  1,  A  to  E).  dl/dV  spectra  measured  with 
the  STM  tip  held  over  the  bare  graphene  near  these  artificial  nuclei  re- 
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vealed  the  emergence  of  a  resonance  in  the  electronic  LDOS  (i.e.,  a  qua- 
si-bound  state)  as  the  number  of  dimers  in  a  cluster  was  increased  from 
one  to  five  (Fig.  1,  A  to  E).  Spectra  acquired  near  2-dimer  clusters  (Fig. 
IB)  displayed  greater  electron-hole  asymmetry  than  single  dimer  spec¬ 
tra,  as  well  as  an  extra  oscillation  in  LDOS  at  high  energies  above  the 
Dirac  point.  Spectra  taken  near  3-dimer  clusters  (Fig.  1C)  exhibited  even 
stronger  electron-hole  asymmetry,  and  the  oscillation  in  LDOS  began 
here  to  coalesce  into  a  resonance-like  structure  near  Fs=  +0.30  V.  For 
four  dimer  clusters,  the  intensity  of  the  resonance  feature  increased 
markedly  and  moved  down  in  energy  to  the  Dirac  point  (Fig.  ID).  For 
the  5 -dimer  cluster,  the  resonance  shifted  below  the  Dirac  point  (Fig. 
IE).  The  formation  of  this  resonance  (or  quasi-bound  state)  as  nuclear 
charge  was  increased  is  the  “smoking  gun”  for  atomic  collapse,  as  we 
will  shortly  describe. 

In  order  to  better  characterize  the  spatial  dependence  of  the  quasi¬ 
bound  state  that  appears  for  high-Z  artificial  nuclei,  we  performed  dlldV 
mapping  at  the  resonance  energy  around  a  5-dimer  cluster  (Fig.  2).  The 
resonance  intensity  displayed  a  highly  symmetric  distribution  around  the 
5 -dimer  cluster  despite  a  marked  asymmetry  in  the  arrangement  of  di¬ 
mers  in  the  cluster.  The  intensity  of  the  resonance  extended  outward 
from  the  cluster  center  to  a  distance  greater  than  10  nm.  Thus,  the  reso¬ 
nance  was  derived  from  charge  carriers  moving  in  the  pristine  graphene 
that  were  spatially  separate  from  the  dimers  making  up  the  artificial 
nucleus. 

The  behavior  of  the  quasi-bound  (atomic  collapse)  state  observed  for 
high-Z  artificial  nuclei  depended  on  whether  it  was  occupied  by  elec¬ 
trons  or  empty,  as  can  be  seen  in  the  gate-dependent  spectra  of  Fig.  3 
which  were  all  measured  3.7  nm  from  the  center  of  a  5 -dimer  cluster. 
Red  arrows  in  Fig.  3  indicate  the  energy  of  the  atomic  collapse  state  at 
each  back-gate  voltage,  while  black  arrows  indicate  the  energy  of  the 
Dirac  point  extracted  by  measuring  dl/dV  spectra  on  graphene  far  away 
(~20  nm)  from  the  cluster  at  each  back-gate  voltage.  In  the  p-doped  re¬ 
gime  (-60  V  <  Vg  <  -15  V)  the  quasi-bound  state  shifted  lower  in  energy 
(relative  to  EF )  as  more  electrons  were  added  to  the  system,  similar  to  the 
shift  of  the  Dirac  point  except  that  the  resonance  transited  from  above  to 
below  the  Dirac  point  as  the  p-doping  was  decreased.  In  the  n-doped 
regime  (0  V  <  Vg  <  +30  V),  the  resonance  shifted  below  EF  and  its  be¬ 
havior  changed  dramatically.  Here  the  resonance  intensity  decreased  and 
essentially  disappeared,  although  a  small  rise  remained  in  the  dl/dV  just 
above  the  Dirac  point. 

The  quasi-bound  state  that  develops  as  the  Ca-cluster  charge  in¬ 
creases  is  the  atomic-collapse  eigenstate  of  a  supercritical  charged  nu¬ 
cleus.  This  explanation  is  supported  by  Dirac  equation-based 
calculations  of  the  expected  LDOS  for  positively  charged  nuclei  with 
different  charge  numbers  (8-10,  14),  which  we  compare  to  our  experi¬ 
mental  spectra.  These  dl/dV  simulations  assume  a  2D  continuum  Dirac 
model  for  undoped  graphene  in  the  presence  of  a  Coulomb  potential.  A 
boundary  condition  at  short  distance  from  the  nucleus  must  be  chosen  to 
account  for  the  finite  size  of  the  clusters.  We  thus  set  V(r)  to  the  constant 
value  V(r0)  for  r  <r0,  a  similar  condition  as  used  in  previous  theoretical 
simulations  of  supercritical  nuclei  in  vacuum  (1,  2).  We  fixed  the  value 
of  r0  to  1  nm,  which  matches  the  geometry  of  our  clusters.  We  found  that 
varying  r0  over  the  range  0.5  nm  <  r0  <  1.5  nm  did  not  markedly  affect 
the  quality  of  the  fits  to  our  data  [simulations  using  other  r0  values  are 
included  in  fig.  S4  (19)].  The  only  essential  fitting  parameter  in  our  sim¬ 
ulation  was  the  ratio  of  the  cluster  charge  to  the  critical  value,  Z/Zc.  Here 
we  define  Z  as  the  “effective  charge”,  i.e.,  the  screened  cluster  charge  as 
seen  by  Dirac  electrons.  The  value  Z  absorbs  the  effects  of  intrinsic 
screening  due  to  graphene  band  polarization  (14)  and  screening  due  to 
the  substrate.  The  critical  value  in  this  case  is  Zc  =  hvF/( 2e2)  ~  0.25, 
where  vF  is  the  Fermi  velocity  [this  Zc  expression  differs  slightly  from 
the  one  used  for  isolated  3D  atoms,  see  supplementary  materials  for  a 
detailed  discussion  (19)].  Additional  perturbations  from  lifetime  broad¬ 


ening  and  electron-phonon  coupling  are  accounted  for  as  in  (14)  and 

(21). 

Our  data  suggests  that  clusters  with  just  one  or  two  Ca  dimers  are  in 
the  subcritical  regime,  and  so  Z/Zc  here  is  determined  by  matching  the 
amplitude  of  the  electron-hole  asymmetry  between  simulation  and  ex¬ 
periment  according  to  the  method  described  in  (14).  Our  clusters  com¬ 
prised  of  3  or  more  dimers  are  either  transitioning  into  (3  dimers)  or 
have  fully  entered  (4  and  5  dimers)  the  supercritical  regime,  as  evi¬ 
denced  by  the  sharpness  and  energy  location  of  the  quasi-bound  state 
resonance  (8,  9).  For  these  clusters,  Z/Zc  is  determined  by  matching  the 
quasi-bound  state  resonance  energy  between  the  simulation  and  experi¬ 
ment.  The  resulting  simulated  dl/dV  spectra  from  the  Dirac  equation  are 
plotted  in  Figs.  1,  F  to  J,  next  to  the  corresponding  experimental  data. 
The  main  features  seen  in  our  experimental  data  are  well  reproduced  by 
the  Dirac  equation  simulations.  In  particular,  our  Dirac  equation  simula¬ 
tions  reproduce  the  increase  in  electron-hole  asymmetry  observed  in  the 
subcritical  regime  as  the  number  of  dimers  increases  from  one  to  two,  as 
well  as  the  emergence  of  the  supercritical  atomic  collapse  state  as  the 
number  of  dimers  is  increased  from  three  to  five.  The  Z/Zc  values  ex¬ 
tracted  from  our  Dirac  equation  fits  are  Z/Zc  =  0.5  ±0.1  (one  dimer),  0.9 
±0.1  (two  dimers),  1.4  ±  0.2  (three  dimers),  1.8  ±  0.2  (four  dimers),  and 
2.2  ±  0.2  (five  dimers)  [further  details  of  the  simulations,  including  bare 
LDOS,  can  be  found  in  the  supplementary  materials]. 

In  order  to  check  that  the  magnitude  of  Z/Zc  extracted  for  Ca  dimers 
from  our  Dirac  equation  fits  is  physically  reasonable,  we  additionally 
performed  a  completely  separate  calculation  of  the  charge  state  expected 
for  a  Ca  dimer  adsorbed  to  graphene  using  ab  initio  density  functional 
theory  (DFT).  This  separate  calculation  (which  had  no  fitting  parame¬ 
ters)  yielded  a  single-dimer  charge  ratio  of  Z/Zc  =  0.6  ±  0.3  [only  the 
single-dimer  case  was  calculated  using  DFT,  see  supplementary  materi¬ 
als  (19)].  This  is  in  agreement  with  the  value  Z/Zc  =  0.5  ±0.1  obtained 
via  our  Dirac  equation  simulations,  and  thus  lends  further  support  to  our 
overall  interpretation  of  the  data. 

We  considered  a  possible  alternative  explanation  for  the  quasi-bound 
state  resonance  seen  for  our  Ca  dimer  clusters  involving  multiple  scatter¬ 
ing  of  electrons  (22)  between  different  Ca  dimers  within  a  cluster.  This 
mechanism,  however,  can  be  excluded  because  the  distance  between  Ca 
dimers  is  d  ~  2  nm,  which  would  result  in  a  peak  energy  measured  from 
the  Dirac  point  of  E  =  hvFk  ~  hvF(n/d)  ~  1  eV.  This  value  exceeds  by 
more  than  an  order  of  magnitude  the  energy  of  the  state  observed  for  Ca 
clusters  on  graphene  (at  most  0.05  eV).  The  large  spatial  extent  (Fig.  2) 
of  the  quasi-bound  state  resonance  shows  that  it  arises  from  an  extended 
Coulomb  field  and  thus  also  rules  out  the  possibility  that  it  is  simply 
caused  by  local  impurity  hybridization  with  underlying  carbon  atoms. 
The  radially  symmetric  distribution  observed  is  consistent  with  the  pre¬ 
dicted  behavior  of  the  atomic  collapse  state  having  lowest  angular  mo¬ 
mentum  (j  =  1/2)  and  lowest  energy  (8,  9).  In  principle,  an  infinite 
number  of  resonances  (corresponding  to  different  principle  quantum 
numbers)  should  appear  in  the  supercritical  regime.  However,  because 
their  energies  follow  an  exponential  sequence  and  their  spatial  extent  is 
inversely  proportional  to  their  energy  measured  from  the  Dirac  point  (8, 
9),  only  the  lowest  energy  state  should  be  detectable  within  our  experi¬ 
mental  resolution  (screening  suppresses  all  states  with  sufficiently  large 
spatial  extent). 

The  strong  doping  dependence  of  the  intensity  of  the  atomic  collapse 
state  (Fig.  3)  can  be  partially  explained  as  a  screening  effect.  In  the  p- 
doped  regime  the  atomic  collapse  state  is  experimentally  observed  to 
shift  to  lower  energy  with  respect  to  the  Dirac  point  for  decreased  p- 
doping.  This  can  be  explained  by  free-electron-like  screening  (23,  24) 
since  reduced  screening  should  result  in  an  increase  of  the  effective 
charge  of  the  artificial  nucleus.  The  n-doped  regime,  however,  shows  a 
completely  different  behavior.  Here  the  amplitude  of  the  atomic  collapse 
state  essentially  disappears  as  it  moves  beneath  EF  (i.e.,  as  it  becomes 
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occupied  with  charge  carriers).  We  believe  that  this  behavior  arises  from 
internal  correlation  among  the  quasi-localized  charge  carriers  that  inhab¬ 
it  the  atomic  collapse  state  (25)  (only  one  out  of  the  four  electron  states 
allowed  due  to  spin  and  valley  degeneracy  can  be  occupied  at  any  given 
time).  Interaction  between  these  electrons  reduces  the  single-particle 
spectral  function  as  measured  in  STM  spectroscopy,  but  the  precise 
mechanism  for  this  reduction  remains  to  be  determined. 
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Normalized  dl/dV 


Fig.  1.  Evolution  of  charged  impurity 
clusters  from  subcritical  to  supercritical 
regime.  (A  to  E)  dl/dV  spectra  measured  at 
different  distances  from  the  center  of  Ca- 
dimer  clusters  (i.e.,  artificial  nuclei) 
composed  of  1-5  dimers  (initial  tunneling 
parameters:  \/s  =  -0.50  V,  /  =  60  pA,  AC 
modulation  Vrms  =  6  mV).  “Center”  here  is 
defined  as  the  average  coordinate  of 
dimers  within  a  cluster.  All  spectra  were 
acquired  at  the  same  back-gate  voltage  (Vg 
=  -30  V)  and  each  was  normalized  by  a 
different  constant  factor  to  account  for 

exponential  changes  in  conductivity  due  to 
location-dependent  tip-height  changes  (14, 
26,  27).  (Insets)  STM  topographs  of 

atomically  fabricated  Ca-dimer  clusters.  (F 
to  J)  Theoretical  normalized  dl/dV  spectra 
(obtained  from  Dirac  equation)  for 

graphene  at  same  distances  from  dimer 
clusters  as  in  (A)  to  (E)  for  the  nuclear 
charges  Z/Zc  of  (F)  0.5  (1 -dimer  cluster), 
(G)  0.9  (2-dimer  cluster),  (H)  1.4  (3-dimer 
cluster),  (I)  1.8  (4-dimer  cluster),  and  (J)  2.2 
(5-dimer  cluster).  Black  dashed  lines 

indicate  Dirac  point,  red  arrows  indicate 
atomic  collapse  state  observed  both  in 
experiment  and  in  Dirac  simulation. 
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High 


Fig.  2.  dl/dV  map  near  a  5  Ca-dimer  cluster  at 
the  energy  of  the  supercritical  quasi-bound  state 
resonance,  as  marked  by  red  arrow  in  Fig.  IE 
(tunneling  parameters:  \/s  =  +0.20  V,  /  =  15  pA, 
Vg  =  -30  V).  The  Ca  dimers  appear  as  slightly 
darker  disks  near  the  center  of  the  dl/dV  map. 
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Fig.  3.  Dependence  of  atomic  collapse  state  on  doping.  Gate 
dependent  spectra  acquired  at  a  lateral  distance  of  3.7  nm  from 
the  center  of  a  five-dimer  Ca  cluster  (initial  tunneling 
parameters:  \/s  =  -0.50  V,  /  =  60  pA,  AC  modulation  \/rms  =  6 
mV).  Curves  are  shifted  vertically  for  clarity.  Red  arrows 
indicate  energy  of  the  atomic  collapse  state  at  each  back-gate 
voltage.  Black  arrows  indicate  energy  of  the  Dirac  point 
extracted  by  measuring  the  dl/dV  spectrum  on  graphene  far 
from  the  cluster  center  (~20  nm)  at  each  back-gate  voltage. 
The  atomic  collapse  state  intensity  is  quenched  in  the  n-doped 
regime. 
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